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Abstract 

We investigate the evolution of the gravitational potential in Rastall scalar 
field theories. In a single component model a consistent perturbation theory, 
Q_j formulated in the Newtonian gauge, is possible only for 7 = 1, which is 

the General Relativity limit. On the other hand, the addition of another 
canonical fiuid component allows to consider the case 7 7^ 1. 
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1. Introduction 

The nature of dark matter and dark energy is one of the most impor- 
tant issues today in physics. There are strong observational evidences in 
astrophysics and cosmology for the existence of these two exotic components 
of the cosmic energy budget, indicating that about 95% of the Universe is 
composed of dark matter (about 25%) and dark energy (about 70%), but no 
direct detection has been reported until now. The usual candidates to dark 
matter (neutralinos and axions, for example) and dark energy (cosmological 
constant, quintessence, etc.) lead to very robust scenarios, but at same time 
they must face theoretical and observational issues. For recent reviews on 
the subject, see for example 

An interesting proposal concerning the nature of dark matter and dark 
energy are the unification models. According to the latter, the whole dark 
sector is a manifestation of a single entity. The paradigm is a perfect fluid 
model called Chaplygin gas but recently it has also been shown that 
viscous models may lead to unification scenarios In spite of their great 
appealing, however, unification models suffer from severe problems when 
confronted with observations, since the parameter estimations from different 
tests lead to contradictory values [sj . One way to surmount this conflict is to 
encode the unification model in a scalar-tensor framework, the exotic fluid 
being described by a self-interacting scalar field. It is not easy to implement 
this idea, since a canonical self-interacting scalar field has a sound speed 
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equal to the speed of light, and it cannot represent dark matter in the past 
evolution of the universe, as required by the unification program j9|. 

A scalar model that is able to represent a realization of dark matter and 
dark energy can be obtained with non-canonical self-interacting scalar fields. 
One example is a scalar field obeying the structure of Rastall's theory, for 



which the usual conservation law for the matter content is modified [10 
The price to pay is the loss of a Lagrangian formulation, at least in the 
context of Riemannian geometry. In Rastall's theory a new dimensionless 
parameter 7 is introduced, which measures the departure from the usual 
equations of General Relativity. When 7 = 1, the General Relativity theory 
is recovered. As in the case of other unification scenarios, the theory is able 
to satisfactorily reproduce the kinematic background observational tests (e.g. 
based on type la supernovae surveys), but essentially reduces to General 
Relativity if a hydrodynamical approach is used for the fluid obeying the 



new conservation laws and the matter power spectrum data are used [11 
However, the agreement improves if a non-canonical scalar field is employed 
instead of a fluid. Moreover, if 7 = 2 this non-canonical scalar field may 



behave as dark matter, and may lead to a unification scenario [12 



Can Rastall's theory pass another important test, represented by the Inte- 
grated Sachs- Wolfe effect? In order to answer this question, the gravitational 
potential $ must be analysed. Using the newtonian gauge, we find an aston- 
ishing property of Rastall's theory: in a scenario with just one component, 
given by the non-canonical Rastall scalar field, only homogeneous solutions 
for $ are admitted, unless 7 = 1. But homogeneous solutions for $ are not 
real perturbations, since they can be re-absorbed in the background metric 
through a suitable reformulation of the background functions. Thus, the 
case 7 = 1 seems to be forced. This fundamental drawback can be cured if a 
two-fluid model is formulated: one scalar field obeying the modified conser- 
vation laws, and a fluid obeying the canonical conservation law. The main 
conclusion is: Rastall non-canonical scalar field requires a matter compo- 
nent (baryons, for example) in order for the theory to make sense at the 
perturbative level. 

In the next section, we introduce the Rastall non-canonical scalar field. 
In section 3, we discuss the sound speed issue in this theory, and in section 
4 we investigate the evolution of the gravitational potential and find the 
constraint 7 = 1. In section 5, we show how to soothe such constraint, by 
adding a canonical fluid component, and we study the gravitational potential 
in a special case. In section 6, we present our conclusions. 
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2. Scalar field in Rastall's theory 

In Rastall's theory |10i] a scalar field (p is characterised by the following 
stress-energy tensor: 

= </',M'^,. - + 9^.(3 - 27)\/(0) , (1) 

where 7 is a parameter. When 7 = 1, we recover the corresponding theory 
in General Relativity. 

We consider a spatially fiat Friedmann-Lemaitre-Robertson- Walker met- 
ric 

ds"^ = dt^ - a^{t)5ijdx'dx^ , (2) 
and its perturbation written in the newtonian gauge 

ds^ = [1 + 2^{t,x%dt^ - a^{t)[l - 2^{t, x')]6,jdx'dx^ , (3) 

where $(t,x*) is the gravitational potential. 

Consider a scalar field perturbation of the form x*) = (f>o{t) + 5(j){t, x^) 
and derive from Eq. ([T]) together with Eqs. ([2]) and ([3]) the background and 
first-order perturbed mixed- component stress-energy tensors. The former is 

(°^n = p = |0o' + (3 - 27)\/(0o) , (4) 
= 0, (5) 
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^0o'-(3-27)F(0o) 



5^- , (6) 



where the dot denotes derivative with respect to t and p and p are the scalar 
field energy density and pressure, respectively. The perturbative quantities 
have the form 



5T% = 5p = 700^0 -7'^>0o' + (3 -27)1^^50, (7) 
5T\ = 4>o6(l),, (8) 

STj = -6p6', = [(7 - 2)00(50 + (2 - 7)$0o' + (3 - 27)^^50] 5', , (9) 

where , i denotes the spatial derivative with respect to the coordinate and 
:= dV {(f)) / d(f) . The modified Klein-Gordon equation has the covariant 
form 

□0+(3-27)K^ = (l-7) '^''/7r > (10) 
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where 00 := 0;^'" . From Eq. (llOp . it appears clearer that 7 = 1 restores the 
General Relativity case. Insert again (p = (pQ + 6(j) into Eq. (ITU!) and employ 
metric ([2]) and ([3]), in order to find 

700 + SH(f>o + (3 - 27) V;^ = , (11) 

where H := a /a is the Hubble parameter. Employing the conformal time 
drj = dt/a{t) we write 

70'o' + (3 - 7)^00 + (3 - 27)aV<^ = , (12) 

where the prime denotes derivation with respect to the conformal time and 
1-L := a' /a. The perturbed modified Klein-Gordon equation has the following 
form: 

7^0" + (3 - 7)^0' - V250 
-(3 + 7)00$' + 2(3 - 27)aV^<l> + (3 - 27)a%<^50 = , (13) 

where V^^ := d'^V {(j)) / d(j)\ 

Differently from the 7 = 1 case, in Rastall's theory we have more degrees 
of freedom. Indeed, considering Tf^^.^ = : 

0,, [00 + (3 - 27) V^] = (1 - 7) , (14) 

one notices that the contraction with 0'^^ gives the Klein-Gordon equation, 
Eq. ( !T0|) . However, Eq. ( !T^ are actually four independent equations. Only 
in the 7 = 1 case they reduce to only one, namely the usual Klein-Gordon 
equation. Using Eq. ( ITU]) into Eq. ([H]) one obtains 

0,.0;a;/30'°0''' ' 0,a0'°0;.;/30''' = . (15) 

At the background level, being 0o dependent only on the time, Eq. (|T5l) is 
identically satisfied. For small perturbations, Eq. ( JTSl) gives the perturbed 
Klein-Gordon equation for z/ = 0, whereas for u = i one has 

(1 - 7)0o50, = (1 - 7)0o50, - (1 - 7)0o'^,i ■ (16) 
For 7 7^ 1 one can cast the above equation as follows: 



= , (17) 

which appears to be an additional constraint that we have to take into ac- 
count together with the Einstein's equations. 
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3. The scalar field speed of sound in Rastall's theory 



In [12| the authors investigate the case corresponding to 7 = 2, which is 
able to reproduce the ACDM scenario both at the background and at the 
perturbative level. A possible explanation for the success of the case 7 = 2 
may reside in the fact that in such instance the speed of sound vanishes, as 
we show now. The speed of sound is defined as the ratio := 5p/5p, which 
is gauge-dependent. Therefore, it makes sense to consider its value in the 
reference frame where the substance whose collapse is being investigated is 
at rest; we denote such quantity as Cg- 

we employ the formula 



Following [13 



5p 



''5p + 3aHp{l + w) (c 



(18) 



which links the pressure perturbations to the energy density ones, both in a 
generic gauge, via Cg^. In this formula, is the adiabatic speed of sound, 
defined as := p/p and which, for the Rastall scalar field that we are 
investigating, has the form 



3g(2- 7)00 + 2(3 -27)\/^ 
3i^70o 



(19) 



where we have employed Eqs. (j4]), and the equation of motion f lTTl) . More- 
over, in formula (|T8|) . w := p/ p, k is the wavenumber coming from a normal 
mode decomposition and 6 is defined via 



a{p + p)9 := d'6T° 



(20) 



Substituting in Eq. flTSj) the expressions for 6p and 6p, that we have found 
in Eqs. ([7]) and (ED, and Eqs. dE]) and ([20]) we obtain 



700^0 -7'^'0o +{3-2-f)V^6(^ + 3H4)o6(f) 



(2-7) 700^0 -T-^-^o +(3-27)V;^(50 + 3i/0o50 



(21) 



which clearly gives that Cg^ = (2 — 7) /7 and therefore the case 7 = 2 implies 
= 0, which is a favouring case for the collapse. 
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4. Evolution of the gravitational potential 



Following , we calculate the Einstein tensor from metric (|3]) and com- 
bine it with the perturbed stress-energy tensor in Eqs. ©-Q. In particu- 
lar, as we implicitly anticipated in writing the perturbed metric ([2]), since 
5T*j oc we have only one gravitational potential. See [13] for more detail. 

We obtain 

v^* - 3n (H$ + $') + 7 (^^ -n')^ = 

4716- [70[,50' + (3 - 27)a%<50] , (22) 
+ = 47rG'0'o<50,, , (23) 

+ sn^' + {2n' + 7^2) $ + (2 - 7) {n^ -h')^ = 

A-kG [(2 - 7)01,50' - (3 - 27)aV<^50] , (24) 

where we have also used the background relation At:G(I)q = TH? — 1-L' . 

Now we reduce the above system in two different ways which, however, 
will give different results unless we choose 7 = 1. Let us employ a normal 
mode decomposition. The second equation can thus be written as 7^$ + $' = 
47rG'0[,(50 and = -k"^. 

Now, if we sum or subtract the first equation with the third, use 7^$+$' = 
47rG0Q50 in order to eliminate 50 and the equation of motion (|T2l) in order 
to eliminate 0o we obtain: 

$" + 3H$' + {2n' + n^)^ = 
^ ~ ^ [-k^^ - m {n^ + $')] - ^^(3 - 27) (n^ + $') , (25) 
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_A;2$ _ 3H^— ^ {m + $') - ^^(3 - 27) {n^ + $') . (26) 

7 700 

These equations can be identical only if 7 = 1, which is the General Relativity 
limit of Rastall's theory. 

Actually, these equations could be consistent also if k = 0, but what does 
this mean? Going back to the Einstein equations, before the normal modes 



7 



decomposition, k = would mean = 0, which imphes that $ should be 
a homogeneous field, i.e. $ = ^{r]). Note that any spatial linear dependence 
for $ would be unacceptable since the field would diverge at infinity. 

But if $ is homogeneous, then we are not facing real perturbations. In 
fact, it is sufficient in metric ([3]) to redefine the time and the scale factor as 
follows: 

dP = [l + 2^t)]dt^ , a\t) = a\t)[l - 2$(t)] , (27) 

and we obtain again the FLRW metric in the usual reference frame, with i 
the cosmic time. 

Therefore we conclude that a consistent perturbation theory, formulated 
in the newtonian gauge, is possible only for 7 = 1. 

Is such result compatible with the new constraint f lT7|) that we have 
showed to exist in Rastall's theory, for 7 7^ 1? The answer is yes, since 
combining Eqs. ( 122|) and (l23|l it is possible to obtain Eq. ( !T7|) . Therefore, 
the latter is not an actual constraint: it is already embedded in Einstein's 
equations. 

5. The role of a fluid component 

The result found in the previous section appears to be different if we 
introduce another component together with the Rastall scalar field. Indeed, 
let us consider a perfect fiuid with equation of state p = wp and w constant. 
We write its total (i.e. background plus perturbed) stress-energy tensor as 
follows: 

T% = p(l + 5), (28) 
= -p(l + w)vi, (29) 
r, = -{p + 6p)6'j, (30) 

where 6 := 6p/p is the usual density contrast, 6p is the density perturba- 
tion, 6p is the pressure one and Vi is the velocity. We assume adiabatic 
perturbations, i.e. 6p = c^Sp, where Cg = w. 

Employing again the normal mode decomposition, we rewrite the system 
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of linearised Einstein equations as follows: 

- - m CH$ + $') + 7 {n^ -n')<^ = AnGa^dp^ + AnGa^dp , (31) 



k (7/$ + $') = A'KGk(t)'Q5ct) + 47rGp(l + w)v , (32) 

$" + 3H$' + 3^2$ - 7 (H^ - ?^') $ = 47rG'a25p^ + 47rGa2(5p , (33) 
where we have defined: 

5p^ := ^ict>'oH' + (3 - 27)V:<^50 , (34) 

Sp^ := ^(2 - 7)0'o50' - (3 - 2^)V^S<P , (35) 

and V is the velocity potential defined by Vi = —v^-Jk. 

With the fluid variables and the relation 5p = c^6p we have a total of four 
unknowns ($, 50, (5p, f ) but only 3 equations. Therefore, it is impossible to 
obtain again a constraint as strong as 7 = 1. 

Let us investigate in some detail the coupled system fluid plus Rastall 
scalar field. Multiplying Eq. (^T\i by 2 — 7, Eq. ([3SD by 7 and subtracting 
the two we obtain: 

7$" + QH^' + - 27 (?^^ - H') $ + (2 - 7)fc^$ = 

-871^0^(3 - 27)y^50 + AnGa^ {^cl + 7 - 2) 5p . (36) 

From this equation we eliminate 5(j) with the help of Eq. f l32|) . obtaining 

7$" + m^' + - 27 (H^ -H')^ + {2- 7)A;2$ = 

+47rGaV (70^ + 7 - 2) 5 , (37) 

where we have used V^^ = V'/(j)'. We now assume the fluid to satisfy its own 
energy-momentum conservation, separately from the scalar field, in order to 
gain one more equation necessary to solve the system. From T'^^-^ = 0, for 
the fluid component only, we get: 

S' = -(I + w){kv - 3^') , (38) 

kr'^ 

v' = -n(l-3w)v + '—6 + k^. (39) 

^ ^ 1+w ^ ' 
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Now we have to specify the background evolution, i.e. the function "H. Its 
general form is 



SttG 



(40) 



In order to simplify considerably Eq. fl37j) . we assume that the potential is a 
constant, i.e. V = 0. Therefore, it is going to play the role of an effective 
cosmological constant. If the potential is a constant, it turns out from the 
Klein-Gordon equation flTT]) that 



0[, = noa-(3-7)/7 ^ 

and the Friedmann equation takes on the following form: 
^2 



-6/7 



with the definitions 



8nG{3-2j)V 



ml 



~3Hr 



(41) 



(42) 



(43) 



Finally, trading the conformal time for the scale factor, we write the coupled 
system of Einstein equations plus the fluid equations as 



n 6 + 7 

7t7 + 

n a 



^6-27 



$ + (2 - 7) 



fc2 



^^?fioa-^(^+"') {icl + 7 - 2) 5 



27/2 

5a = -(I + W) 



kv 
T-La 



(44) 
(45) 



1 kr"^ k 

va = —{1 - 3w)v + —-^6 + —<!> , (46) 
a [l + w)Ha Ha 

where the subscript a denotes derivation with respect to the scale factor. 
We start the evolution from a, = 10~^ and choose as initial conditions 



0, 



With the choice w = and 7 = 2, Eqs. ( 142|) and ( l44ll reproduce the same 
dynamics of the ACDM model. It is curious that this seems to happen only 
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Figure 1: Evolution of $ and S for the choice w = and flo = 0.04. The black solid lines 
represent the case 7 = 2, i.e. the ACDM one. The red dashed lines correspond to 7 = 1.9 
whereas the blue dot-dashed ones to 7 = 2.05. We have chosen a scale k = h Mpc~^. 




Figure 2: Evolution of $ and 5 for the choice w = and VLq = 0.04. The black solid lines 
represent the case 7 = 2, i.e. the ACDM one. The red dashed lines correspond to 7 = 1.5 
whereas the blue dot-dashed ones to 7 = 1. We have chosen a scale k = 10^^ h Mpc^^. 
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when a standard fluid is added to the Rastall scalar field. In Figure [T] and [2] 
we plot the evolution of $ and 6 for the choice w = and Qq = 0.04. That 
is, we are assuming that the perfect fiuid under consideration is a baryon 
component. We have also chosen a scale k = 10~^ h Mpc~^. It seems that, 
for 7 > 2, the growth of the density contrast of the fiuid component is 
enhanced. This is probably due to the fact that Cg becomes negative and 
therefore the collapse of the scalar field is unimpeded. For 7 < 2 the growth 
of S is sensibly hampered. It is curious in Figure |2] how the gravitational 
potential suffers a larger decrease for the case 7 = 1.5 than for the 7 = 1 
one. We would have expect the contrary, since for 7 = 1 the speed of sound 
is Cg = 1, whereas for 7 = 3/2 the speed of sound is Cg = 1/3, i.e. smaller. 
Note that such discrepancy is not present in the plots for 6, i.e. the growth 
for 7 = 1.5 is larger than the one for 7 = 1. Therefore, such effect is probably 
due to the different background evolution. 

For completeness, we display here also the evolution of 6^ := 5p^/p^, 
which can be easily calculated from Eq. fl3T]) . once we know $ and 5. The 
results are plotted in Figure |3i 




0.0 0.2 0.4 0.6 0.8 1.0 



Figure 3: Evolution of for the choice w = Q and JIq = 0.04. From top to bottom: 
7 = 2.05, 2, 1.9, 1.5, 1. We have chosen a scale k = 10"^ h Mpc"^ 

Some comments about Eq. ( |T7|) are in order, since the latter establishes a 
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strong connection between perturbations in the scalar field and in the gravita- 
tional potential. For a single scalar field component, Rastall's theory reduces 
to General Relativity, thus everything runs as in the standard lore. On the 
other hand, the coupling with matter brings new features. Here also comes 
the complexity given by the new conservation law of Rastall's theory, which 
admits many consistent alternatives. If the matter component conserves sep- 
arately (as we have investigated in this paper), Eq. flT7|) remains untouched. 
However, there are examples in which it may change, e.g. the case in which 



matter exchanges energy with the scalar field or the one investigated in [15 
where two fiuid components were considered, one of them still conserving 
separately, whereas the conservation of other depending on the curvature. In 
this case, it is not difficult to show that relation ( 1T7|) gives place to 



4>o J '2 



rendering the situation somewhat different. In [15[ other possible couplings 
between the two components are also evoked, which may lead to other vari- 
ants for Eq. ( |T7|) . 

Even if our interest in this work is the late-times universe, we may ask for 
potential consequences of Eq. ( 1T7|) and its possible variants for the primordial 
spectrum. One point is that Rastall's scalar model requires another compo- 
nent in order to make sense. But, an adiabatic primordial spectrum can be 
naturally implemented mainly if the matter component is subdominant with 
respect to the scalar component. The isocurvature component can also be 
implemented, in principle, since it requires a zero total (scalar plus matter) 
density fluctuation 5ptot = 0, and this even if the relation f|T7|) still holds, as 
in the case where the other component conserves separately, without direct 
interaction with the scalar fleld. When other types of interactions between 
both components are considered, as in the equation above, the isocurvature 
perturbation can still exist. But, in general, the detailed predictions for the 
spectrum must differ from the standard cases, mainly in the isocurvature 
case. This may open a new path of investigation concerning the speciflc 
predictions of Rastall's theory for the primordial spectrum of perturbations. 

Equation ( IT71) also reminds the relation between the spatial curvature and 



the inflaton fleld, in the standard inflationary scenario ^6[. However, such 
connection still has to be investigated in detail, in order to understand how 
inflation could be implemented into Rastall's theory. 
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6. Conclusions 

The scalar formulation of Rastall's theory of gravity may allow a consis- 
tent unification of dark matter and dark energy for the background evolution 
of the universe. We have shown in this paper that, on the other hand, its 
single component version (with the non-canonical Rastall scalar field as the 
only matter content) is perturbatively inconsistent: the compatibility of the 
perturbed equations requires a homogeneous gravitational potential, which 
is equivalent to a redefinition of the background functions and not to real 
perturbations. The Rastall scalar theory may admit consistent perturbative 
scenario if another (canonical) fluid component is added. 

For a two-fluid model, we numerically evaluate the behaviour of the grav- 
itational potential and that of the density contrast for the scalar field com- 
ponent. For some cases, as 7 ~ 1.5, a behaviour very similar to that ob- 
tained in the General Relativity case with a quintessence field is obtained. 
Although this does not represent an exhaustive study of the Rastall two- 
component model, the results here reported indicate that consistent scenar- 
ios may emerge from a Rastall unification model for dark energy. We hope to 
present in future a more exhaustive study with a detailed comparison with 
cosmological observational data. 
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